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Abstract 

The Wheeler-DeWitt equation of Friedmann models with a massless scalar quan- 
tum field is formulated with arbitrary factor ordering of the Hamiltonian constraint 
operator. A scalar product of wave functions is constructed, giving rise to a prob- 
ability interpretation and making comparison with the classical solution possible. 
In general the behaviour of the wave function of the model depends on a critical 
energy of the matter field, which, in turn, depends on the chosen factor ordering. 
By certain choices of ordering the critical energy can be pushed down to zero. 

1 Introduction 

We consider homogenous, isotropic cosmological models with a scalar, massless 
field as source. Both the gravitational degree of freedom and the scalar field are 
quantized, thus leading to a simple quantum mechanical system. The notoriously 
ambiguous factor ordering in the Hamiltonian constraint operator in the Wheeler- 
DeWitt equation is left open by introducing a two-parameter family of factor or- 
derings. In all the cases we arrive at a Bessel or a modified Bessel equation as wave 
equation for negative and positive curvature, respectively. A probability interpre- 
tation of the wave functions is supported by considerations about the possibility 
to derive an equation of continuity from the Wheeler-DeWitt equation, which is 
infiuenced by the factor ordering. Then for not too small values of the scale factor 
of the cosmological model the probability densities, derived in a natural way from 
the wave functions, approach their respective classical limits. 

An analysis of the factor ordering issue in closed Friedmann models can be 
found in That work is aimed at a critical review of the Hartle-Hawking no- 
boundary PI and the Vilenkin tunneling wave functions ^ of the universe. In 
the meantime progress has been made in Loop quantum cosmology based on a 
quantum theory of gravity which allows to calculate and make physical sense 
of wave functions at and beyond the classical cosmologic singularities and where 
the factor ordering plays a key role, too. Whereas in pP the potential of the scalar 
field is approximately constant and its kinetic energy does not play a role in the 
Wheeler-DeWitt equation, we assume a free, massless field, only coupled minimally 
to gravity, leading to separable, exactly solvable equations. By introducing two real 
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parameters (and in one case complex ones) we consider a larger degree of freedom 
of factor ordering than where one integer parameter is used. 

In all the cases, i.e. for the closed, the open, and the spatially flat model, there is 
a critical energy of the matter field with the associated gravitational wave function 
being regular (or going to infinity) near the singularity for energies below it and os- 
cillating infinitely many times for energies above it. The existence of a critical energy 
of the matter field is also observed in the framework of Loop quantum cosmology 
Here, in the simpler formalism of standard quantum cosmology, where we have 
reduced the system to a mechanical one before quantizing, but introduced factor 
ordering parameters, we can control the factor ordering by varying the parametrs 
arbitrarily. By variation of the factor ordering it is possible to shift the critical 
energy, it is also possible to choose such orderings, where the critical energy goes 
down to zero. So, if one would prefer a unique behaviour of the wave functions at 
the singularity at all energies, e.g. to satisfy a DeWitt condition of zero probability 
at the singularity, this would restrict the possibilities of factor ordering. 



2 The classical model 



We start with the action 



c I -Lg + Lm 1 dn 

K 



where Lm is the lagrangian for a scalar field 



Lm 



(1) 



(2) 



and Lg is the geometrical lagrangian. Greek indices are spacetime, V{(f)) is the 
potential of the scalar field. Under the assumption of homogeneity and isotropy 
a standard simplification of variables leads to the Robertson- Walker metric (RW 
metric) and after integration of the action ((2j over a unit volume^, we get the 
hamiltonian by a straightforward Legendre transformation. 



H = N 
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+ ka 



(3) 



^The Ricci scalar R is only a function of the time for the RW metric, so only y/—g is a function of 
spatial variables. In the closed case we can integrate over the whole volume, in the open case we must 
integrate over some finite, spatially constant "test volume" , so we can study time dependence of the scale 
for this volume. After this step, the lagrangian is only a function of time and the equations of motion 
deduced from this lagrangian by variation are equivalent with the Einstein equations for the RW metric. 
With these assumptions we have finite dimensional mechanical lagrangian and hamiltonian theory for 
the RW models and we can use this in quantum theory. In this example the reduction to a mechanical 
system before quantization works, but not generally 10;^ . 
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Here A'^ is the lapse function, a the scale parameter, c the speed of light, k Ein- 
stein's gravitational constant k = and k is the normalized sectional curvature, 
which can be +1,-1 and for the closed, open and flat case. Momenta are defined 
in the standard way 

dL a?^ 

The set of Einstein's equations following from the action ^ is nonlinear, which 
is a good reason for using the hamiltonian formalism for the solution of this problem. 
In the massless case the equations are 



1 fPaV , 3 



= -24 If J ■ w 
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" ~ 12 a ' 

= 0, (8) 



(7) 
(8) 
(9) 



Together with the hamiltonian constraint 

H = (10) 
we have a system equivalent with the Einstein equations for the RW metric. 

3 Classical solutions 
3.1 Solution for k = 1 

The momentum of the massless scalar field is constant. Thus one can write condition 
()1U() in the form 



p, = ±12a^/(^)'-l, (11) 

where a. is the maximal radius of the univerl, given by = g and obtained for 
d = 0. If we combine condition ([TT|l with the dynamical equations (6) and (7), we 
get 



Or, 



a2 ^ (13) 
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Pa = 12 +1 , (12) 
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The expansion speed a in the second equation is the first integral of the first equation 
(if we put Pa = —12aa). The second equation leads to the integral 

± / dt, (14) 



where we suppose < a < a^- The solution can be written in form of a GauB 
hyper geometric function [S] 



/ x3 /lQ'7 / \4N 



For the inversion of this function, leading to a(t), we have to consider the ex- 
panding and the contractiong branch seperately. Denote the time of the maximal 
radius as tm = Fc{a.m)- The expanding branch is obtained by choosing the positive 
sign on the 1. h. side of ([T5|) and to = 0, whereas the negative sign and to = 2trn. 
yields the contracting branch. After the unification of the two branches we get the 
typical time evolution of the closed cosmological model without cosmological con- 
stant shown in figure 1, i. e. the time evolution from the big bang to the collapse. 
For completeness, the time of the maximal radius is in SI units 



r(|)r(|) 
r(f) 



C^ra = ^,5' ~ 0, 599 am- (16) 



Our exact solution is an explicit confirmation a theorem in [H] for the case of the 
massless scalar field. 

On figure n the evolutions of the scale a for the closed and the open (bold line) 
model, with Om = 1, p = 1 (see definition below eq. ((T^ l and c = 1 are displayed. 

Note that (|15j) is a solution of a GauB equation which has second linearly inde- 
pendent solution: a constant. The solution for the field is 



o ; / \ 4 



,o±y-arctanhyi-^-j , (17) 
so we have fully described our system on the classical level. 

3.2 Solutions for k = —1 and = 

The open case is analogous. The only difference is a plus sign under the square 
root in 1)11(1 and in the integral ()14() and of course there is no maximal radius of the 
universe. The solution is 

±(*-^o) = f (^) '2^1(^4- i'-(^) l-^o(a), (18) 
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Figure 1: The time evolution of the RW models. 
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where = is the formal analog of am- If we choose = for the positive sign 
in H18|) we get a solution for the scale parameter increasing in time. In this solution 
the expansion speed of the universe is converging in time to a constant, this follows 
from the equation 

.^ = (f)Vl. (20) 

Because a is increasing forever, it follows from this equation that d ~ 1 for 
large t. This is typical for open cosmological models with vanishing cosmological 
constant. 

In the flat case A; = we obtain the simple solution 



a{t) = const. Vi. 



(21) 
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4 Wheeler-DeWitt equation 



One of the possibilities how to choose a representation by selfadjoint operators on 
£^{R+,^;R,d(j)) is [H] 

ail){a,<j)) := a'il^{a,(j)) TTTp(a,cl)) := — i^a ^^^"' '^^ (22) 

oa 

(j)'4){a, 4>) := (l)'4){a, 4>) P4>ip{a, (/)) := -i/i^^^^-^. (23) 

For the motivation for this (noncanonical) choice of vr and the above Hilbert space, 
in which it is selfadjoint, see fSJ and quotations therein. 

The Poisson brackets lead to the affine commutation relations 

{a,7r} = a — > [a,7r]=iM, (24) 
{4>,P4>} = [4>,P4>] =^fi, (25) 

where classically tt = apa- 

With the usual substitution pa —ihd/da in the hamiltonian ^ we obtain the 
following two-parameter family of Wheeler-DeWitt equations^ 



6-ka + ca'^V{4>)] 'il^{a,<l)) =0 (26) 



24 c a* da daa'^ 4 c d(f>'^ k 

with the real numbers i,j, k fulfilling i+j + k = 1. This is our approach to the factor 
ordering problem. It is easy to show that the factor ordering problem is impossible 
to resolve by means of the semiclassical limit, because by assuming ifj = exp j^S 
one obtains a Hamilton-Jacobi equation without regard to the numbers i,j-,k after 
neglecting terms containing h. 

There are some orderings which one can prefer. Most popular is to consider 
expressions with derivatives as a Laplace-Beltrami (LB) operator 1^ in superspace. 
There is only one such ordering: i = 2,j = —1, A: = ie. ■^'^'Ti'- This is denoted as 
"D'Alembert operator ordering" in ^j. Another choice, namely i = k = 0, j = 1, 
was used by Vilenkin jHI- 

On the other hand, there exists an infinity of possibilities of ordering. We will 
show solutions for a massless field in all these cases. Nevertheless, someone might 
prefer special selections: for example all symmetric orderings of a and vr are 

^ l^d_^±^_ .27) 



^AU quantum formulae are in SI units. 
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where n — 1 = j, or in terms of the canonical operators a and pa 
p2 I d I d I 



(28) 



This operator is not selfadjoint on the Hilbert space mentioned before, but on 
/^^(M"*", da; M, d(/)). A complex generalization (here i is the imaginary unit) of a 
special case of (^7)) is 

TT^ 1 / 1 _ _ 1 \ 1 . +.i . 1 

TTTT + TTvr^ = — — -7ra=^'^'7r^ -. (29) 



a3 2 a? J af(i±i) a^{i±^) 

Also this last choice represents manifestly selfadjoint operators on £'^(M"'", ^). 

For the massless scalar field the WD equation ()26() is simply solvable for every 
ordering and curvature, because one obtains an equation with separable variables. 
An Ansatz ip{a,(j)) = A{a)ip{(j)) leads to 

1 d 1 d 1 , 144 , A2 , , 

r- T:A{a) - - — -TrTT^kaAia) + -^Aia) = 0, (30) 

a' da ai daa^ ^ ' {STrlj^ ^ ^ ^ ^ ' ^ ' 
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2 +—^{^) = 0, (31) 



here Ip is the Planck length. The solution of the second equation is a plane wave (p ~ 
g±ip^0^ which corresponds to the result of the classical theory, where the momentum 
of the scalar field is an integral of motion. 

Now we introduce a new variable b = (^f'^ and a new function B{b) as 
A{a) = a^^ ~ B{b{a)), so equation (|Hn|) takes the form 



B' f~ v^^ 
B" + --(k + ^]B = (32) 



with 



' - [(j + 1)^ - • (33) 

V may be real or imaginary in dependence on the relation between the energy or 
momentum of the field (j), given by A = \J^P(f), and the ordering parameter j. In 
all the cases we can write the a-dependent part of the wave function as 

with By being a solution of H32() . 
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5 Inner product 



If we introduce the notation = (x^ = a^x^ = (p), we will get after some calcula- 
tions from (|2()j) 



5^ (l^* 



with the flow vector 



■= r - ^ 9^ J ' ^ 5^ - ^ 50 J • 

The metric in x" space is Gap = diag(— 24a, 4a'^), which follows directly from 
If we choose i = k, then equation 1)35^ represents an equation of continuity for the 
flow Ja- Of course, the right hand term of (|35|) and the first term of (|36|) vanish, 
if the o-dependent part of the wave function -0 is a real function of a real variable. 
This holds for the closed model, as we shall see in the next section. 

In the symmetric case i = k and in the case of a real wave function of a a volume 
integral of (|35|) in two-dimensional superspace spanned by a and <j) leads to an 
integral over the flow through some closed loop (= one-dimensional closed surface) 
in (two-dimensional) superspace, which is equal to zero. Moreover, in the closed 
case there is only one nonzero component of J"^, namely J^, flowing in the direction 
of (j). Thus, if (p is considered as internal time, for the closed model the integral over 
a is conserved in time and represents an indefinite norm corresponding to the inner 
product = J ^{ip*d^ip — (fd^ip*). (As long as we have a certain value of p^, 

or the parameter A, respectively, the role of the internal time <p is analogous to the 
role of time in an energy eigenstate in quantum mechanics - there is no dynamical 
evolution. Dynamics arises when we take superpositions into account.) 

The sign of is arbitrary because the conservation law holds for the positive 
as well as for the negative sign of p^ (if i = k) and, as long as p^ ^ and as there 
are no superpositions of wave functions with positive and negative values of Ptf,, we 
can interprete \'ip*d^'ip — ipd^jyip*] as probability density, which will be seen to be in 
accordance with the classical one. 

The requirement of writing (|26() in covariant form, i. e. with the derivatives in 
form of the LB operator ^vrvr mentioned above, leads to the formula 



V„ J" = ^= daiV^-n = 0, (37) 




with the same as in (|36jl and with an interpretation as probability density for 
\/ —GJa rather than for the components Ja- 

Now we can integrate the equation (|37() over a volume in superspace and using 
the Gaufi integral theorem, which holds in all dimensions, we get 



Gd2x= / (9„(^/^J")d2x = / ^/^Jan'^dl, (38) 
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where is the hne element of the closed loop and is the normal to this loop. 
The conserved quantity is \/ —GJa- The DeWitt condition A(a 0) is postulated 
for a function which has an interpretation as a probability density. This condition 
means that the classical singularity at a = is avoided. In our case it means, 
that |J2(a 0)1 = for the symmetrical ordering or \\/—GJ2{a ^ 0)| =0 for 
the LB ordering. But we also must study whether the right hand side integral 
is correctly defined. The overlap of these two conditions also restricts the factor 
ordering parameters. 



6 Closed model 



6.1 Solutions of the hamiltonian constraint 

If we choose k = 1, we get from (|32|) a modified Bessel equation of the order v, 
which can be viewed as a Bessel equation with purely imaginary argument. The 
general solution has the form 

B{h) = cih{b)+C2K^{b). (39) 

Here Iy{b) = i^'^Ju{ib). This function is not limited in the asymptotic domain and 
for this reason it is not normalizable. The second function is a Basset function 
which is defined as 

K^ib) = ^iei-/7i(i6). 

Let's first assume that the order v is imaginary and write i/ = le 
case when the square of A is above a critial value, coming from (|33j) . 

In terms of the absolute value of the canonical momentum this means 

>Pcr = y¥|U + l|. (42) 

If we put b = x^, we can rewrite the Basset function of imaginary order in the 
form 



(40) 
This is the 



(41) 



TT 



sinhTre V 2 

m=0 



(3 cos ( e In 



a sm 



e In ■ 



(43) 



where a = a(m, e) is the real and (3 = /3{m, e) is the imaginary part of r(l+m)r(l + 
m + I/). From this it follows that the Basset function of pure imaginary order and 
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with arguments from is a real function of a real variable. In the asymptotic 
limit Basset functions of real or imaginary order reduce to the form 

^.W^^yje-'- (44) 

On the other hand the function is limited for — > 0+, as can be seen from (|43|) . 
For these reasons we will accept the Basset function of imaginary order as the 
appropriate part of (|39|l . So the physically acceptable solution of the Wheeler- 
DeWitt equation (|Hn|l in the closed case, for the momentum of the scalar field above 
the critical value, is 

A,(a) = c(A)(y/^^)'+'?i^i,(„(^(^)^). (45) 

6.2 Inner product 

Taking the scalar product in the Hilbert space 

£2(M+, d^;M,d(/.) from the fore going 
section, we arrive, up to a A-dependent factor, at the following norm of a mode 
^A(a) 



2+k-i 1^2 

ie\ 4tt\ 

This is finite if2 + A; — 'i>0. If2 + A; — i>l, then the integrand goes to zero for 
a ^ 0. Interpreted as probability density to find the model at a certain value of a, 
the DeWitt condition of the avoidance of the classical singularity a = 0, is satisfied. 
The above conditions depend on the factor ordering parameters i and k; for the 
selfadjoint choice H27() . where k — i = 1, they are satisfied. 

The alternative choice (|28j) with i = k and the corresponding integration mea- 
sure da in the Hilbert space leads to exactly the same result. For both choices 
normalizability restricts factor ordering, the DeWitt condition enhances the restric- 
tion. Figure 2 compares the assumed quantum probability density with the classical 
one, obtained from (fT^ as the inverse of the velocity d, 

Pc^ = l = -r^,- (47) 

In the case of undercritical field momentum, when < (^^^ , we obtain a 



Basset function of real order, K^, ( ^ ( f" ) ) ' 'with 



-V(i + 1)2-4A2 (48) 



as a-dependent part of the wave function. Obviously the domain of A of these 
solutions depends on the factor ordering parameter j, if j = —1, there is only 
one solution Kq for A = 0. Given the case that j ^ —1 there are solutions for 
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A G {—^^^, ^-^), but they may be normalizable or not. From the behaviom' of the 
Basset functions of order u ^ for a close to zero, 

^'^fcfelV""^ (49) 



it follows that the normalizing integral A\ (a) exists if 



l + k-i - yj{3 + IY - 4A2 > -1. (50) 

For the LB operator, for example, with i = 2, j = —1, and k = 0, the only solution 
of this type with = is not normalizable. Here the normalizable solutions with 
imaginary u for > and those with < are separated by one non-normalizable 
solution with = 0. 

A normalizability condition resulting from (|5()j) is 

i-k<2. (51) 

If we consider factor ordering parameters close to this boundary, characterized by a 
small number e, i — k = 2 — e, we can express i and k by j and a further parameter 
6, using i + j + k = 1, 

i = l{3-j-6), k = ^{-l-j + 6). (52) 

Now we insert back into ()5U|) and obtain a lower boundary for A, 



2 / r\ 2 

< A^ (53) 

So, for overcritical values of A^, there are normalizable, oscillating wave func- 
tions, for undercritical values there may be, depending on the parameter 6 = 2—i+k, 
normalizable non-oscillating functions for 

Whether or not there are normalizable wave functions for undercritical values of 
A, depends on the relation between the deviance of i — k from 2 and j. The same 
applies for the ocurrence of non-normalizable functions, which seperate normalizable 
functions with positive and negative values of A from each other. Of course, for the 
symmetrical factor orderings the functions are normalizable. 

6.3 The DeWitt condition 

In terms of the field momentum we can separate solutions, which fulfil some conti- 
nuity equation, into two parts. The first kind of states, for which holds 



'^iO-+i)V4^i(^-+i))' (55) 
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will be called dynamical (overcritical) . The remaining states will be called non- 
dynamical (undercritical) . We called these states dynamical, because close to the 
singularity the probability densities derived from these states are comparable with 
the classical ones. Moreover, dynamical states must be represented by a function 
oscillating near the classical singularity. Effectively the classical quantity 1/d ~ 
is comparable with the corresponding quantum expression near the classical singu- 
larity (for solutions fulfilling some continuity equation), as one can see from (j43|) 
and H45() . Unfortunately, the mathematically nice Laplace-Beltrami factor ordering 
has no nondynamical states for non zero field momentum (see (|55|) with j = — 1). 
Any other factor ordering has nondynamical states for ^ 0. 

Every dynamical state, which fulfils a continuity equation, fulfils the DeWitt 
condition. Nondynamical states fulfil the DeWitt condition iff 



p<i, e 



If j = 1 the unification of dynamical and nondynamical states satisfies the 
DeWitt condition for momenta G M\{0}. Note that the symmetric ordering 
j = 1 is derived as semiclassical limit in loop quantum cosmology 9 . It has a 
maximal set of values of the field momentum, for which J2(a — > 0) = holds. 
Any other choice fulfilling a continuity equation (i.e. symmetrical) doesn't fulfil the 
DeWitt condition in some finite momentum interval. 
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Figure 2: Classical and quantum probability density. 
On figure |2l classical (grey) and quantum (black) probability densities for j = 1 
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are displayed. We can see a density constructed for a dynamical (bold lines; = 
2\/6) and a nondynamical state (thin lines; = 1/3). 

On the classical level there is no reason for a restriction of A. This could be 
an argument in favour of a factor ordering with j = — 1, n = or j = =p3i — 
1. Specifically this holds for the LB operator ^tttt. On the other hand, in loop 
quantum cosmology ^ other factor orderings are preferred and an analogous critical 
field energy (of the order of the Planck energy) for a dynamical evolution of the 
model exists. 

Nevertheless, an argument in favour of the symmetrical l|28|) or LB ordering is 
that for it we can correctly formulate the continuity equation. LB ordering may be 
interpreted as ordering of the selfadjoint operators ^tttt. 



7 Solution for the open model 



In spite of principal problems with the canonical quantization of open cosmological 
models |1U| , here we simply take over the Wheeler-DeWitt equation from the closed 
model and change the sign of the curvature and justify this step by the correct 
classical limit of the resulting wave function. For k = —1 (jSH) becomes a Bessel 
equation and as in the foregoing case we must distinguish between the cases z/^ > 
and < 0. For v = ie in the case > A^j. we get a solution in terms of two 
linearly independent Bessel functions of imaginary order, 



B{b)=CiJUb) + C2J-ieib), 

as before. (The two Bessel functions are complex conjugate, 
solution can be given in terms of the Hankel functions Hj^^^ and H^"^^ 
function for the open model becomes 



(56) 

Alternatively the 
The wave 



Axia) 



l+(k-i)/2 



Cl(A)Jie(.)U(^) +C2(A)J-ie(.,^(t 



(57) 



For convenient assumptions about the factor ordering, i.e. k — i = for the 
Hilbert space with the measure da, and k — i = 1 for the measure da/a, the square 
of the absolute value of the above functions approximates very closely the classical 
probability density l/a = a^ j a^ + p^ from (|2()|). see fig. 4. 

For real u the solution of H32|) may be written as superposition of Hankel func- 
tions Hi^\b) and H^\b) or the corresponding Bessel functions of the first and 
second kind. In the first case only the expression b\Hu\b)\'^ as probability density 
would be constant for large a and so display the classical behaviour, but it diverges 
at a = for v > 1/2. If we admit only J^(6) for > we have a function that 
vanishes at the classical singularity, but it does not have the character of an out- 
going or ingoing wave function for large a. Moreover, the absolute squares of the 
Bessel functions of real order are oscillating and do not approximate the classical 
behaviour shown in fig. 4. 
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Figure 3: Plot of | Ji,(a2)|2. 

In the spatially flat case, when k = 0, the situation is analogous. For i/ = ie 
the solutions of (|32|) are ~ a^^'^ and have a dynamical interpretation as in- and 
outgoing wavefunctions. Given an appropriate factor ordering with respect to the 
chosen Hilbert space measure, they are normalizable to a (^-function. For real v the 
solutions do not have these properties. 

8 Summary 

In the spatially flat case there is a simple distinction betwee over- und undercritical 
values of the field momentum p^p, or the parameter A, respectively: For too small 
values there is no dynamical wave function, corresponding to the classical solution. 
In the other cases there may be a nonvanishing range of the above parameters be- 
tween the critical value and zero, in which the wave functiones might be physically 
acceptable. For k = —1 this means there are two linearly independent solutions 
iif(^) and H^"^^ and a probability density, according to the classical limit. In the 
closed case A; = 1 we mean by acceptable wave functions normalizable ones. The ex- 
istence and the extension of a range of undercritical energies with normalizable wave 
functions depends crucially on the factor ordering and is more restricted by classi- 
cal considerations in the open case, where the limit of the probability density for 
large a singles out symmetric orderings. A physical interpretation of "undercritical" 
normalizable wave functions is still missing. 

The flow Jq has two nonvanishing components, but we reserve an interpretation 
as probability density for \VGJ2\ in the case of LB ordering, or for the component 
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I J2 1 in the case of symmetric ordering. 

To compare once more with []P, not the form of the potential, but the kinetic 
energy of the scalar field plays a role in our paper. Due to this, rather than "tunnel- 
ing" and "oscillating" regions in superspace, i.e. a different behaviour of the wave 
function in different domains of the scale parameter, we obtain different regions 
concerning the energy of the field modes. Our result is not a rigid restriction of 
factor orderings or the picking-out of some distinguished ones, but rather a subtle 
dependence of the critical values of the field energy on the factor ordering. 

Acknowledgement: Supported by the Ministry of Education of the Czech Repub- 
hc under the project MSM 0021622409 
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